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Abstract 

We considered N x N Wishart ensembles in the class Wq C^n , M) (complex 
Wishart matrices with M degrees of freedom and covariance matrix Sjy) such that 
A^o eigenvalues of S^r is 1 and A^i = — A^'o of them are a. We studied the limit as 
M, N, No and A^i all go to infinity such that ^ ^ c, ^ ^ /3 and < c, /3 < 1. In 
this case, the limiting eigenvalue density can either be supported on 1 or 2 disjoint 
intervals in M_|_. In the previous paper ^28j, we studied the universality in the case 
when the limiting eigenvalue density is supported on 2 intervals and in this paper, we 
continue the analysis and study the case when the support consists of a single interval. 
By using Riemann-Hilbert analysis, we have shown that under proper rescaling of 
the eigenvalues, the limiting correlation kernel is given by the sine kernel and the 
Airy kernel in the bulk and the edge of the spectrum respectively. As a consequence, 
the behavior of the largest eigenvalue in this model is described by the Tracy- Widom 
distribution. 



1 Introduction 

Let X be an M X (assuming M > N) matrix with i.i.d. complex Gaussian entries whose 
real and imaginary parts have variance | and zero mean. Let S^r be an A^ x A^ positive 
definite Hermitian matrix with eigenvalues ai,...,aj\r (not necessarily distinct). In this 
paper, we will consider the case where S^v has only 2 distinct eigenvalues, 1 and a such 
that A^^i of its eigenvalues are a and N — Ni of them are 1. We will assume that ^ ^ c and 
^ — > /3 as A^, M — > cxD and that < c, /5 < 1. To be precise, we will assume the following 

cM - N = n = 0{l), N(3 - Ni = T2 = 0{l), M,N,Ni^oo. (1.1) 

1 1 
Let be any Hermitian square root of Sat. Then the columns of the matrix XE^ are 

random vectors with variances ^y/Oj- Let the matrix he the following 

i?^ = ^4^tx4, (1.2) 
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Then Bn is the sample covariance matrix of the columns of XS^, while Sat is the covariance 
matrix. In particular, Bm is a complex Wishart matrix in the class Wc {^n, M). 

The sample covariance matrix is a fundamental tool in the studies of multivariate 
statistics and its distribution is already known to Wishart at around 1928 (See e.g. PSj) 

V{Bj,) = le-^^f-(^-^^^) (det S)""'-^ , M > N, (1.3) 

for some normalization constant C. 

Let yi, . . . ,yN > he the eigenvalues of the sample covariance matrix B]\f. Then in 
the case where A^^i eigenvalues of S^v is a and N — Ni is 1, the joint probability density 
function (j.p.d.f) for the eigenvalues of Bjy is given by 

N 

ny) = -Y- n(^^ - n yf'"" \yT"^''"'^^'"\,^. _ , (1-4) 



where Zm,n is a normalization constant and , aj are given by 

df = j, a, = 1, l<j<iV-iVi, 

df = j-N + Nu aj=a, N - Ni < j < N. 

In this paper we will study the asymptotic limit of the Wishart distribution with ^ — ^ c 
and ^ — /? as M, N ^ oo in such a way that < /?, c < 1. In this case, the empirical 
distribution function (e.d.f) F/v of the eigenvalues will converge weakly to a nonrandom 
p.d.f. F, which will be supported on either 1 or 2 intervals in IR+. By applying the results 
of [IS] to our case, we can extract properties of the measure F from the solutions of an 
algebraic equation (See Section [3] for details) 

za^^ + {A2Z + ^2)^' + {z + 5i)e + 1 = 0, 
A2 = {l + a), 52 = a(l-c), (1.5) 
fii = 1 -c(l -/?) +a(l -c/?). 

The results in [15] imply that the real zeros of the function ^^^^p- determines the boundary 

points of the support of F. Since the zeros of coincide with the zeros of the following 
quartic polynomial, 

a'{l - c)e + 2{a\l - c/3) + a(l - c(l - 

+ (1 - c(l + a\l - cP) + 4a)e' + 2(1 + a)^ + 1 = 0, ^ ' ' 

the real roots of (11. 6p are important in the determination of Supp(F). In particular, we 
have the following result. (See Theorem E]) 

Theorem 1. Let A be the discriminant of the quartic polynomial ( fi.6]) . If A < 0, then 
the support of F consists of a single intervals. 
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In our previous paper [2H|, we have shown that Supp(-F) consists of 2 disjoint intervals 
if and only if A > and have obtained the asymptotic eigenvalue statistics in that case. 
In this paper, we consider the case when A < and together with [28], we have proven the 
universality in this class of complex Wishart ensemble for all A 7^ 0. When A = 0, a phase 
transition takes place and the support of the eigenvalues splits into 2 disjoint intervals. 

When A < 0, we also have the following expression for the density function of F (See 
Theorem [6]) . 

Theorem 2. Let A he the discriminant of the quartic polynomial 1^1. 6{) . Suppose A < 
and let 71 < 72 be the 2 real roots to U.6\) and 73, 74 be the 2 complex roots. Let X^, 
k = 1, . . . ,4 be the following 



1 1-/3 
— + c- 

7k 1 + 7fc 



+ c 



a/3 



1 + a7fc • 



k = l 4. 



Then all the Xk are distinct and the p.d.f F is supported on [Ai,A2] with the following 
density dF{z) = p{z)dz 



3 

2^ 



r[z) + 



27a4 



r{z) 



1 



;i-7) 



where D^{z) and r{z) are given by 



D,{z) = {l-afX{{z-X,\ 
r[z) = — I —z + 



+ 



27 \ a3 

9^2 2Al 



6A2BI ^ _2 



z-^ + 



9B2 9B1A2 27 6AIB2 



-1 



The constants Ax, Bi and B2 in the above equation are defined by 

A2 = (l + a), B2 = a{l-c), Bi = l-c{l- (3) + a{l-c(3). 

The cube root in ^3. 23\) is chosen such that \I~A G M for A G M and the square root is 
chosen such that \I~A > for A > 0. 

Remark 1. Since D^lz) can be written as 



-27{riz)Y-A{p{z)r 



'11 
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for some polynomial p{z) in z"^ . We see that ifr{z) vanishes at any of the Xk, then D^{z) 
will have at least a double root at these points, hence r(Afc) 7^ 0. From this and ( [i. we 
see that the density p{z) vanishes like a square root at the points Xk, k = 1,2. 

p{z) = P^\z - Xk\-^ + O {{z - Xk)) , z-^Xk. (1.9) 

TT 

An open problem in the studies of Wishart ensembles is the universality and the dis- 
tribution of the largest eigenvalue. Although the Wishart distribution is known for a long 
time, results in the universality and the largest eigenvalue distribution were only obtained 
recently and only for spiked models [21] whose covariance matrices are finite perturbations 
of the identity matrix [6], [7|, [8], [IE], [2lj, [22j, [21], [25], [32], [33]. The resuh in this 
paper is one of the few results obtained for models with covariance matrices that is not a 
finite perturbation of the identity matrix. (See also [2D], in which the largest eigenvalue 
distribution was also derived for ensembles whose covariance matrix is not a finite pertur- 
bation of the identity matrix. However, in [2Q], the parameters c = /5 = ^ and the 



M' N 

eigenvalues of the covariance matrix have to satisfy a condition which is not true in our 
case.) 

In this paper, we use an important result by Baik, Ben-Arous and Peche [7] which shows 
that the correlation functions of the eigenvalues can be expressed in terms of a Fredholm 
determinant with kernel KM,N{x,y). In [H] and [16], the authors have expressed this 
kernel in terms of multiple orthogonal polynomials (See Section [2] for details) and have 
shown that the m-point correlation function for the Wishart distribution (11.31) is given by 

7^(f '^)(i/i, ■■■.y^) = det {KMAy^.yk)),<,,k<m (i-io) 

where 7lm^'^\yi, . . . ,ym) is the m-point correlation function 

7^L*^'^)(2/l,••■,yJ = 7T7^^ / V{y)dy„,+,...dy^. (1.11) 

(iV-m)!JK+ Jr+ 

By computing the asymptotics of the correlation kernel in (ll.lOp through the asymptotics of 
multiple Laguerre polynomials, we have proven the universality of the correlation function 
when M, N and Ni ^ 00. 

Theorem 3. Suppose A in TheoremU\is negative. Let p{z) be the density function of F 
in Theorem\^ Then for any xq G (Ai, A2) and m & N, we have 



M^oo\Mp{xo)J " V Mp(xo)' Mp{xo)J 



N,M^oo \Mp{xo 

( sin 7r(uj — Uj 



det 



'1.121 



uniformly for any {ui, . . . , Um) in compact subsets of. 
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On the other hand, let xq = Xk, k = 1,2, then for any m E N, we have 



' Ai{ui)Ai'{uj) - Ai{ui)Ai{uj 



m 

I A1I7/..1A1 I7/„1 — Al I 7/..- 1 All 7/.„- 1 \ 

det 



'U'i III 'j / ■ ■ 1 



;i.i3) 



uniformly for any {ui, . . . , Um) in compact subsets o/M™, where Ai{z) is the Airy function 
and pk, k = 1,2 are the constants in ( (i.gj) . 

Recall that the Airy function is the unique solution to the differential equation v" = zv 
that has the following asymptotic behavior as z — > cxo in the sector — vr + e < aig{z) < vr — e, 
for any e > 0. 

1 2 3 / _3 \ 

Ai{z) = re 1 + 0(z 2) j ^ _7r + e < arg(z) < tt - e, 2:^00. (1-14) 

2^zi ^ ^ 

where the branch cut of z^ in the above is chosen to be the negative real axis. 

Since the limiting kernel takes the form of the Airy kernel (I1.13p . by a well-known result 
of Tracy and Widom [HI] , we have the following 



Theorem 4. Let 7/1 he the largest eigenvalue of B^, then we have 

lim P Uyi - A2) {MP2Y <s)= TW{s), (1.15) 

where TW{s) is the Tracy-Widom distribution 

TW{s) = exp (t - s)q'^{t)dt^ , (1.16) 

and q{s) is the solution of Painleve II equation 

g"(s) = sg(s) + 2g3(s), 

with the following asymptotic behavior as s ^ 00. 

q{s) ~ — Ai(s), s — >■ +00. 

The results obtained in this paper are obtained through the Riemann-Hilbert analysis. 
As in [12] and [26], a Riemann surface of the form (11.51) . together with the zero set of a 
real function h{x) related to this Riemann surface are essential to the implementation of 
the Riemann-Hilbert analysis. This function h{x) is given as follows. If we express the 
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Figure 1: Possible sheet structures of the Riemann surface. For the implementation of 
the Riemann-Hilbert analysis, we need to show that the Riemann surface has the sheet 
structure shown on the right hand side. 

solutions ^ in (11. 5p as analytic functions of 2;, then the three solutions to (11.51) behave as 
follows when z 00. 

^i{z) = --^ + 0{z-^), z-^00, 

6W = -l + ^^^^ + 0(z-^), z^oo, (1.17) 

U^) = -- + ^ + 0{z-'), z^oo. 
a z 

Then the function h{x) is defined by 

h{x) = Re (^^ ^z) - U^)dz^ (1.18) 

In order to implement the Riemann-Hilbert analysis, we must determine the sheet structure 
of the Riemann surface (11.51) and the topology of the zero set of h{x). Since our model 
depends on three parameters c, a and (3, while the models in |T2j and [2B] depend only 
on one parameter a, the determination of both the sheet structure of the Riemann surface 
and the topology of ^ are considerably more difficult in our case and a large part of this 
paper is devoted to resolve these difficulties so that Riemann-Hilbert analysis like those in 
[12] and [26] can be applied. 

To determine the sheet structure of (11.51) . note that although we assume the discriminant 
A in Theorem [T] is negative, which means that the Riemann surface (II. 5p has 2 real and 2 
complex branch points, it is unclear which Riemann sheet do these branch points belong 
to and these different situations will result in different sheet structures of the Riemann 
surface as indicated in Figure [H In order to determine the sheet structure of the Riemann 
surface, we need to analyze the analyticity of the solutions in the vicinity of the points 
Afc in Theorem [21 This requires the type of analysis used in [27] , which is very difficult to 
carry out in our case. As in [28], we overcome these difficulties by showing that the Stieltjes 
transform of the limiting eigenvalue distribution satisfies (11.51) . Then by using properties 
of the Stieltjes transform obtained in [3J, [4J, [5], [I5], [30], we were able to determine the 
sheet structure of the Riemann surface (II. 5p . 
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The determination of the topology of the set in this case is also much more compli- 
cated and a thorough analysis of this set making use of properties of harmonic functions 
is carried out in Section I3.2.1I 

For theoretical reasons, the model studied in this paper is crucial in understanding of 
the phase transition behavior that occurs in Wishart ensembles. (See [7]). When the 2 
intervals in the support of F closes up, a phase transition takes place and the correlation 
kernel will demonstrate a different behavior at the point where the support closes up. With 
the Riemann-Hilbert technique used in this paper, such behavior can be studied rigorously 
as in pTSj and the eigenvalue correlation function near the critical point will be given by 
the Pearcey kernel. For practical reasons, many covariance matrices appearing in fields 
of science are not finite perturbations of the identity matrix. In fact, covariance matrices 
that have groups of distinct eigenvalues are accepted models in various areas such as the 
correlation of genes in microarray analysis and the correlation of the returns of stocks in 
finance. 
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for the financial support provided by the grant EP/D505534/1. 

2 Multiple Laguerre polynomials and the correlation 
kernel 

The main tool in our analysis involves the use of multiple orthogonal polynomials and the 
Riemann-Hilbert problem associated with them. In this section we shall recall the results 
in [H] and [16] and express the correlation kernel Km,n{x, y) in (11.101) in terms of multiple 
Laguerre polynomials. In Section HI we will apply Riemann-Hilbert analysis to obtain the 
asymptotics of these multiple Laguerre polynomials and use them to prove Theorem [SI 

We shall not define the multiple Laguerre polynomials in the most general setting, but 
only define the ones that are relevant to our set up. Readers who are interested in the 
theory of multiple orthogonal polynomials can consult the papers [T], [2], [E], [23] • Let 
-^"1,^2(2^) be the monic polynomial such that 




(ix = 0, i = 0, . . . 



X 



dx = 0, i = 0, . . . , n2 — 1. 



,ni 



-1, 



(2.1) 



and let Qm ,112(2^) be a function of the form 




(x)e-^- + <^,„^(x)e 



(2.2) 
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where (^) (^) polynomials of degrees ni — 1 and n^ — \ respectively, 

and that QuY^n^i?^^ satisfies the following 

^ Qni,n2{^)^ dx = ^ 2 = m + n2 — 1 (^•^) 

The polynomial L„^^„2(^) is called the multiple Laguerre polynomial of type II and the 
polynomials Al^_^ (x) and [x) are called multiple Laguerre polynomials of type I 

(with respect to the weights x^'^~^ e~^'^ and x'^'^~^ e~^''"' [1], [2]. These polynomials 
exist and are unique. Moreover, they admit integral expressions [Ti] . 
Let us define the constants hn},n2 

and hni,n2 

to be 

L„,„,(x)x-+*^-^e-*^^rfa;, 

^ni,n2 ~ / -^ni,ri2 (■^)"'' ^ ^ 
^0 



Then the following result in and [16] allows us to express the correlation kernel in 
(11.101) in terms of a finite sum of the multiple Laguerre polynomials. 

Proposition 1. The correlation kernel in KM,N{x,y) I^1.10\) can be expressed in terms of 
multiple Laguerre polynomials as follows 

N-AI 

{xy) 2 [x - y)KM,N{x,y) = Lno,Ni{x)Qno,Ni{x) 

"'No,Ni J t \r\ /■ \ 

- -(Y) Lno-i,nA^Nno+i,nA^) . . 



— LNo,Ni-l{x)QNf^^Ni+l{x) 



"'No,Ni-l 



where No = N - Ni. 

This result allows us to compute the limiting kernel once we obtain the asymptotics for 
the multiple Laguerre polynomials. 



3 Stieltjes transform of the eigenvalue distribution 

In order to study the asymptotics of the correlation kernel, we would need to know the 
asymptotic eigenvalue distribution of the Wishart ensemble (11.31) . Let Fn{x) be the empir- 
ical distribution function (e.d.f) of the eigenvalues of (II. 2p . The asymptotic properties 
of -F/v(x) can be found by studying its Stieltjes transform. 

The Stieltjes transform of a probability distribution function (p.d.f) G{x) is defined by 

poo 

mciz) = / T dG{x), zeC-^ = {zeC: lm{z) > 0} . (3.1) 
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Given the Stieljes transform, the p.d.f can be found by the inversion formula 

1 r'' 

G([a,6]) = -hm / Im {mds + ie)) ds . (3.2) 

The properties of the Stieltjes transform of -F/v(x) has been studied by Bai, Silverstein and 
Choi in the papers [15], [3], [5], [30]. We will now summarize the results that we need 
from these papers. 

First let us denote the e.d.f of the eigenvalues of S^r by Hjy, that is, we have 



Furthermore, we assume that as N ^ oo, the distribution Hn weakly converges to a dis- 
tribution function H. Then as — >■ oo, the e.d.f Fn{x) converges weakly to a nonrandom 
p.d.f F, and that the Stieltjes transform mp of F{x) satisfies the following equation [15] . 



Let us now consider the closely related matrix 

= ^XS^Xt. (3.4) 

The matrix B_j^ has the same eigenvalues as -Bat together with an addition M — N zero 
eigenvalues. Therefore the e.d.f Fjy of the eigenvalues of are related to F/v by 

X 

Zat = (1 - Cn)I[Q,oo) + CatFat, Cat = — . (3.5) 

where /[o,oo) is the step function that is on R_ and 1 on IR+. In particular, as X — > oo, 
the distribution Fjy converges weakly to a p.d.f F_ that is related to F by 

F = (1 - c)/[o,oo) + cF (3.6) 
and their Stieltjes transforms are related by 

1 - c 

mF_{z) = V c'mF{z). (3.7) 

Then from (13.31) . we see that the Stieltjes transform rriF^z) satisfies the following equation 



This equation has an inverse [1], [5] 



z{mp) 



1 



+ c 



/ ^ dH{t) 

1 + tmp 



(3.9) 



The points where 



are of significant interest to us as they are potential end points 



of the support of F_, due to the following result by Choi and Silverstein. 
Lemma 1. [ij5j If z ^ supp(F), then m = mF_{z) satisfies the following. 

1. m G M \ {0}; 

2. -^i supp(if); 

3. Let z he defined by ^3. 9^ . then z'{m) > 0, where the prime denotes the derivative with 
respect to mF_ in ^3.9\) . 

Conversely, if m satisfies 1-3, then z = z{m) ^ supp(F). 

This lemma allows us to identify the complement of supp(F) by studying the real points 
m such that z'{m) > 0. 

Remark 2. As pointed out in f^, LemmaUl applies to any distribution G whose Stieltjes 
transform mciz) Ii3. 1\) satisfies an equation of the form 



for some constant cq and distribution Hcit). 

3.1 Riemann surface and the Stieltjes transform 

We will now restrict ourselves to the case when the matrix S^r has 2 distinct eigenvalues 
only. Without lost of generality, we will assume that one of these values is 1 and the other 
one is a > 0. Let < /5 < 1, we will assume that as — oo, A^i of the eigenvalues take the 
value a and Nq = N ~ Ni oi the eigenvalues are 1 and that ^ (3. That is, as oo, 
the e.d.f Hn{x) converges to the following 



By substituting this back into (13. 9p . we see that the Stieltjes transform ^{z) = mp^z) is a 
solution to the following algebraic equation 




dHN{x) dH{x) = (1 - p)6i + 135a. 



(3.10) 




1 



l-(3 



+ o 



a(5 



(3.11) 



1 + e 
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Rearranging the terms, we see that ^ = mp^ solves the following 

zai"" + {A2Z + ^2)^' + {z + B{)i + 1 = 0, 
A2 = (l + a), 52 = a(l-c), (3.12) 
Bi = l-c(l-/3)+a(l-c/3). 

This defines a Riemann surface £ as a 3-folded cover of the complex plane. 

By solving the cubic equation fl3.12p or by analyzing the asymptotic behavior of the 
equation as 2; — > cxd, we see that these solutions have the behavior given by fll.l7p as 
z 00. On the other hand, as z — > 0, the 3 branches of C,{z) behave as follows 

Uz) = -^ + 0{l), z^O, 

^piz)=R, + Oiz), z^O, (3.13) 

e^(z) =i?2 + 0(z), z^O. 

where the order of the indices a, (3 and 7 does not necessarily coincide with the ones in 
(11.171) (i.e. we do not necessarily have a = 1, (5 = 2 and 7 = 3). The constants Ri and R2 
are the two roots of the quadratic equation 

a(l - c)x^ + (1 - c(l - /3) + a(l - c(3))x + 1 = 0. (3.14) 

The functions ^j{z) will not be analytic at the branch points of £, and they will be dis- 
continuous across the branch cuts joining these branch points. Moreover, from (13.131) . one 
of the functions ij{z) will have a simple pole at z = 0. Apart from these singularities, 
however, the functions S,j{z) are analytic. 

3.2 Sheet structure of the Riemann surface 

As explained in the introduction, the determination of the sheet structure of the Riemann 
surface (11.51) involves difficult analysis of the behavior of the ij{z) at the points A^. In 
order to implement Riemann-Hilbert analysis to obtain the asymptotics of the kernel in 
(11.101) . we must show that ii{z) is not analytic at the real branch points Afc, k = 1, 2. This 
will be achieved by making use of the properties of the Stieltjes transform. 

In this section we will study the sheet structure of the Riemann surface C As we shall 
see, the branch ^i{z) turns out to be the Stieljes transform mp_{z) and its branch cut will 
be the support of F. 

By Lemma [U the real points at which ^ > characterize the end points of Supp(F). 
The determine the support, let us differentiate (13.111) to obtain an expression of ^ in terms 
ofe 



dz 



a^(l - c)i^ + 2{a^{l - c(3) + a(l - c(l - /5))^^ 



di e(i + o'(i + <)' 

+ (1 - c(l - /3) + a^{l - cp) + 4a)e' + 2(1 + a)^ + 1 



(3.15) 
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In particular, the values of ,^ at ^ = correspond to the roots of the quartic equation 



a\l - c)e + 2(a^(l - c/3) + a(l - c(l - m .3 

+ (1 - c(l - /3) + a2(l - c/3) + 4a)e' + 2(1 + a)^ + 1 = ^ ' ^ 



Let A be the discriminant of this quartic polynomial, then when A < 0, the equation 
(13.161) has 2 distinct real roots 71 < 72 and 2 complex roots 73 and 74 = 73. One can check 
that the coefficients of (13.161) are all positive and hence 71 < 72 < 0. 
Let Afc be the corresponding points in the 2;-plane 

X, = -- + c^^ + c-^, fc=l,...,4. (3.17) 

7fc 1 + 7fc 1 + «7fc 

Then we have the following 

Lemma 2. The points Xk, k = 1, . . . ,4 are all distinct. 

Proof. The derivative has simple zeros at the points {Xk,7k), k = 1, ...,4, on the 
Riemann surface defined by (13.121) . Therefore, as a function of z, 2 branches of the function 
C,{z) behaves as 

^(z) = 7fc ± Cfc(z - Afe)2 + 0(z - Afc), z ^ Xk, 

near the points Xk- Let i j, then Since 7^ 7^ jj, if Aj = Xj, there will be 4 distinct 
solutions ^{z) to the equation (I3.12p in a neighborhood of the point Aj = Xj, which is not 
possible. Therefore Aj and Aj are distinct. □ 

In particular, the points A3 and A4 are complex. 

Lemma 3. The points A3 and A4 are not real. 

Proof. Since 73 = 74, we see that A3 = A4. If A3 and A4 are real, then we will have A3 = A4. 
This contradicts Lemma [2] and hence the points A3 and A4 are not real. □ 

We will label the points 73 and 74 such that Im(A3) > 0. 

Note that, from the behavior of z{^) in (13. lip , we see that near the points —1 and — -, 
the function z{^) behaves as 

^(0 = ^^^ + 0(1), ^--1, 

. (3.18) 

^(0-^ + 0(1), 

1 + at, a 

and hence z'{^) is negative near these points. From this and (13.150 . we see that z'{C,) > 
on the intervals (—00,71), (72,0) and (0, 00) and none of the points —1 or — ^ belongs to 
these intervals. On (71,72), z'{^) is negative. 

The images of these intervals under the map (13. lip then give us the complement of 
supp(F) in the 2;-plane. Let us study these images 
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Lemma 4. The intervals (—00,71), (72,0) and {0,oo) are mapped by z{C,) to (0,Ai), 
{X2,oo) and (— oo,0) respectively. 



Proof. Since none of the points —1, — ^ and belongs to these intervals both z(^) and z'{^) 
are continuous on these intervals. Moreover, z{^) is strictly increasing on these intervals. 
Therefore the images of these intervals are given by 

2; ((-00, 71)) = (2;(-oo),z(7i)) = (0, Ai) 
^((72,0)) = (2;(72),2;(0-)) = (A2,oo) 
^((0,00)) = (^(0+),^(oo)) = (-00,0). 

where the ± superscripts in the above indicates that the function is evaluated at ±e for 
e ^ 0. □ 

Therefore the complement of supp(F) is given by (recall that F_ has a point mass at 0) 

supp(F)" = (-00, 0) U (0, Ai) U (A2, 00). (3.19) 

Let us now show that Ai < A2. This would imply the support of F is non-empty and 
consists of one interval. 



Lemma 5. Let Ai and A2 be the points defined by 1^3.17^ , then Ai < A2 and hence the 
support of F_ consists of a single interval. 

Proof. Suppose Ai > A2. Let Zq G (A2,Ai). Since both the points —1 and — ^ belongs 
to (71,72), the function z{^) is continuous in (—0x0,71) and (72,0). By Lemma HI these 
two intervals are mapped onto (0, Ai) and (A2,cxd) by the map z{^). Hence there is at 
least one point in each of the intervals (—00,71) and (72,0) that is being mapped onto 
Zq. On the other hand, let Si and S2 be the points in {—1, —\} that is closer to 71 and 
72 respectively, then z{^) is continuous in (71, Si) and (52,72)- Since these intervals are 
mapped onto (— oo,Ai,) and (A2,oo) respectively by z{^), there is at least one point in 
each of these intervals that is mapped onto zq. This would result in 4 distinct points being 
mapped onto zq by the map z{^)^ which is not possible. Therefore we have Ai < A2. □ 

Therefore we have the following 

Theorem 5. Let A be the discriminant of the quartic polynomial 

a\l - c)e + 2(a^(l - cP) + a(l - c(l - P))e 

+ (1 - c(l ~(3) + a^{l - 0(3) + 4a)e' + 2(1 + a)^ + 1 = 0. ^ ' ' 

then if A < 0, the support of F_ consists of a single interval. 
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We can treat (13.121) as a polynomial in ^ then A^, k = 1, ... ,4 are the zeroes of its 
discriminant D^lz) = {az)^ nj<7 (^* " 0)^- 



D3{z) = (1 - ayz* + {2AlBi + 2^252 - AA^ - 12a5i + 180^2)^^ 

+ {Bl + AlBl + 4^251^2 - 12^^52 - 1205^2 + 180^2 + 180^2^1 - 270^)^=^ 
+ {2BiBl + 2A2B2BI - I2A2BI - ABfa + 18aBiB2)z + B^B^ - AB'l 

(3.21) 



The zeros of (I3.2ip then correspond to the branch points of the Riemann surface C. These 
branch points are given on C by (Afc, 7^), for = 1, . . . , 4. 

Since the leading coefficient of D-j,{z) is (1 — a)^ > 0, we see that the sign of D^^i^z) and 
hence the 3 roots of the cubic polynomial (I3.12p behave as follows for z G M. 



1. 2; G M \ [Ai, A2], D^^z) > ^ ^ has 3 distinct real roots 

2. 2; G (Ai, A2), D^i^z) < ^ ^ has 1 real and 2 complex roots. 

In particular, since the roots coincide at the branch points, the 7^, k = 1, . 
values of the double root of the cubic (13.121) when z = Xk- 
We can now compute the probability density F. 



(3.22) 
, 4 is the 



Theorem 6. Suppose A < 0. Then the p-d.f F_ is supported on [Ai, A2] with the following 
density dF{z) = p{z)dz 



3 

2^ 



riz) + 



1 



D.(z) 



^ r{z) 



1 



D.(z) 



\ 



(3.23) 



where D^lz) is given by 1^3. 2 1\) and r{z) is given by 



riz 




2Bl n 

-^z + 



2Ar 



9BIB2 6A2B 



z-' + 



9B2 9B1A2 27 6AIB2 
a? a? a 



The cube root in ^3. 23\) is chosen such that \I~A G M for A G M and the square root is 
chosen such that \I~A > for A > 0. 



The proof of this theorem is the same as Theorem 6 in [28] . 

Let us now show that the function C,i{z) is in fact the Stieltjes transform mp_{z). From 
the asymptotic behavior of the ^j{z) (I1.17P and the fact that the Stieltjes transform mp_{z) 
is the unique solution of (I3.12p that vanishes as 2; 00, we see that 



(3.24) 
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For c < 1, it was shown in that F has a continuous density and hence the Stieltjes 
transform mp^z) does not have any poles. Therefore, by (13. 7p . we see that mp^z), and 
hence ii{z), has the following singularity at z = 0. 

ei(^) = -^ + 0(l), z^O. (3.25) 

Since mp_{z) is the Stieltjes transform of a measure supported on the real axis, it is analytic 
away from the real axis and hence by (13.241) . ^i{z) is analytic at the points A3 and A4 with 
a branch cut on [Ai, A2]. Since A3, A4 are not branch points of the function ^1(2;), they must 
be branch points of the functions ^2 (-2) and ^3(2;). This determines the branch structure of 
the Riemann surface C The branch cut of ^2(2;) and ^3(2;) will eventually be chosen to be 
a contour that goes from A4 to A3 intersecting the real axis at a point in (Ai, A2), but in 
the next section it will be chosen in a few different ways according to the situation. 



3.2.1 Geometry of the problem 

As explained in the introduction, the determination of the zero set of h{x) in (ll.lSp is 
considerably more difficult than it is in [12] and [26]. In this section we will carry out a 
thorough analysis of this set and determine its topology. 

For the implementation of the Riemann-Hilbert analysis, it is more convenient to con- 
sider a measure F/v instead of the measure F_. 

Let Cat = ^ and H]^{t) be the e.d.f 

dHn{t) = {l-f3N)8i+(3NSa. 

where (3^ = jf- 

Then Fjy is the measure whose Stieltjes transform mjv(z) = mp^{z) is the unique 
solution of 

z{mN) = — + CN [ —\ dHN{t), (3.26) 

in C+ that behaves like — ^ as z — 00. Note that, as pointed out in [8], the measure F/v is 
not the eigenvalue distribution for finite N , instead, it is only defined through the equation 
(13.261) . From (13.261) . we see that mis[{z) is the solution of the algebraic equation 

zae + {A2Z + B^)e + {z + B^)i + 1 = 0, 
B^ = a{l-CN), 5f = l-c;v(l-/3^)+a(l-Civ^7v). 

that behaves as — ^ as z ^ 00. If we denote by if {z) the solutions of (I3.27P with asymptotic 
behavior (I1.17p . but with c and {3 replaced by cn and /?iv, then we have {z) = rriM^z). 
In [15], it was shown that the measure has a continuous density on and a point 
mass of size 1 — cat at and hence {z) also have the asymptotic behavior (I3.25P near 
2; = with c replaced by Cat- By Remark [21 all the results in the previous section will 
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remain valid for fl3.27p and with P and c replaced by Pn and cat. We will denote 

the Riemann surface defined by f l3.27p Cn- 

Let Aiv be the determinant of the quartic polynomial fl3.16p 



a\l - cn)^' + 2{a\l - cnPn) + a(l - c^(l - 

+ (1 - cn{1 - Pn) + a'(l - cnPn) + 4a)e' + 2(1 + a)^ + 1 = ^ ' ' 



Since the determinant A of the quartic polynomial (13.161) is continuous in the parameters 
a, P and c, for large enough N, M and A^i, we can assume that the determinant An < 0. 
Then let < 7^ be the real roots of fl3.28p and 7^^ and 7^ be the complex conjugate 
roots of (13.280 and let \^ < and A^, A^ be their images under the map z(^) given in 
(I326D. 

For the time being, we will choose the branch cut of ^i{z) to be the interval [X^ , A^], 
and the branch cut between A^ and A^ to be a simple contour C that is symmetric with 
respect to the real axis and oriented upwards. Across C, the two branches ^2^(2;) and ^^{z) 
change into each other. The branch cut C will be chosen such that it intersects the real 
axis at exactly one point x* that is not equal to or A^. 

We will now define the functions Oj^{z) to be the the integrals of ^f{z). 

Oii^) = r ^ii^)dx, e^{z) = r ^^ix)dx, 9^{z) = r ^^ix)dx, 

Jx^ (3.29) 

ifx*<A^, 1 = 2; ifx*>X^, 1 = 1. 

The integration paths of the above integrals are chosen as follows. If x* < A^, then the 
integration path will not intersect the set CU(— 00, A^) and if x* > A^, then the integration 
path will not intersect the set C U (Af^, cxd). 

Then from (11.170 . (I3.13P and (13.250 . we see that the integrals (I3.29P have the following 
behavior at 2 = 00 and z = 0. 

e^iz) = -log z + + 0{z-'), z^oo, 
<(z) = -(1 - cn) log^ + O (1) , z-^ 0, 
e^{z) = -z + cn{1 - Pn) \ogz + l^ + [z-^] , z^oo, 

e^iz) = oii), z^o, 

e^iz) = -^ + cnPn log z + l^ + (z-^) , ^ ^ 00, 

e^iz) = oii), z^o. 



(3.30) 



for some constants , and . 
Then the set defined by 

Sj = {zeC\ Re{e^{z) -e^{z)) =0} (3.31) 

is important to the Riemann-Hilbert analysis. Let us now study its properties. First note 
that it is symmetric across the real axis. 
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Lemma 6. Let the branch cut C between \^ and be a simple contour that is symmetric 
with respect to the real axis. Then the set S) in l\3. 31\) is symmetric with respect to the real 
axis. 

Proof. Let us consider the transformation z and S,{z) y-^ S,{z) in f l3.27p . This gives 

zalijf + {A2-Z + B^)l{zf + {z + B^)l{-z) + 1 = 0. (3.32) 

This means that if ^{z) is a solution to f l3.27p . then so is iil). Since the functions {z) is 
analytic away from CU [A{^, A^], it must equal either of the functions {z), ^2 i.^) i.^)- 
By considering the behavior of (^) near z = 00 using f ll.l7p . we see that (^) = if{z) 
for i = 1, 2, 3. Let Zq G S), then we have 



Re 



J\^2i^)-C^{^))dx^ =0. (3.33) 



By taking the complex conjugation of 03.331) and making use of the fact that (z) = ^fi^) 
for j = 2, 3, we obtain 

Re (^£i^2i^) - ^3i^))dx I = 0. (3.34) 

Let us now show that Re (^2^(A^) — ^3^(A^)) = 0. Consider an integration contour F from 
A^ to A^ that is symmetric with respect to M. Then we have 



\N \ / „\JV 



Re I / ^ iC^ix) - i^{x)) dx]=Re\ / ' {^^ {x) - (x)) dx 

-Re| / ■ {^^(x)-^^ix))dx 



N I \ I \N 

3 / \ "^4 



where we have used A^ = A4 in the above. Hence we have Re (6*^ (A^) — 0^{X^)) = 0. 
This implies the following 

Re - ^3i^))dx] = Re (£i^2i^) - ^3i^))dx^ = 0- (3-35) 



Therefore if 2:0 G io, then zq is also in Sj and hence Sj is symmetric with respect to the real 
axis. □ 

We will now show that the set Sj is independent on the choice of the branch cut C. 

Lemma 7. Let the branch cut C be symmetric with respect to the real axis. Then the set 
is independent on the choice of the branch cut C. 
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Proof. Let us consider the boundary values of 62 {z) and 0^{z) along the branch cut C. 
Let ^2±i^) ^3±i^) be the boundary values of ^^(2;) and ^^{z) on the left and right 
hand sides of C, then we have C,2±{z) = ^3zf{z) for z E C. Therefore, for z in the upper 
half plane, we have the following 

where the integration is performed along C. Therefore in the upper half plane, the function 
62 (z) — 6^ (z) changes sign across the branch cut C. Hence in the upper half plane, the 
zero set of Re (z) — 9^ (z)) is independent on the choice of C. Now by Lemma the 
set is symmetric with respect to the real axis for any choice of symmetric branch cut 
C. Therefore the set in the lower half plane is just the reflection of in the upper half 
plane which must also be independent on the choice of C. □ 

We will now show that the real parts of ■^2^(2;) and C,^{z) will coincide exactly once on 
the real axis. 

Lemma 8. Let x* be the intersection between C andM, then the real function Re (^^(-z) — 
is continuous on {—oo,x*) and (x*, cxd) and it vanishes exactly once at a point l G M\{2;*}. 

Proof. The location of x* is immaterial as the function Re (^^(-z) — {z)) only changes 
sign across the point x* and hence its zeros on R are independent on the location of x* . 
For definiteness, let us assume that X* <\^. 

Let /i = 2, /2 = 3 for a > 1 and Zi = 3, /2 = 2 for a < 1. From the behavior f lL17p 
of if {z) near z = ±00, we see that there exists R > such that C,f{±R) > (±-R) > 
C,l^{±R). Note that il^{z) and (z) are continuous on (— oo,x*) U (A^, 00). Since these 
intervals do not contain any branch point of fl3.27p and both ^j^ {z) and (z) are real on 
them. Re (^/^ (z)) and Re (z)) cannot coincide on these intervals. In particular, the 
order ^[^{z) > il^{z) must be preserved in (— oo,x*) U (A^,oo). Since if{z) and ii^{z) 
interchange across the branch cut C, we must have if{z) > (z) in the interval (x*, Xf). 
This means that if^{z) > if^{z) on the left hand side of [Af,A^], while if^{z) < il^{z) 
on the right hand side of [\f, A^] and hence Re {^fi^)) -^^ (^3^ (-2)) niust coincide at 
least once in [Af , A^]. We will show that they can only coincide once within [Xf , A^]. 

Inside [Af^,A^], the function ii{z) and one other root if{z) becomes complex and 
are conjugate to each other, while another root ^^{z) remains real. Note that, since 
Cf±{z) = (,fzf{z) on the branch cut [A^,A^], we see that the real part of C,f{z) and 
Cf{z) is continuous across [Xf, A^], while ^^{z) has no jump discontinuity across [Xf, A^]. 
Since neither if{z) or ^^{z) are equal to the branch ii{z) in {Xf , A^), the real function 
R-^ {if{z) — iniz)) must equal either of ±Re (^^(-z) — ^^iz)) and hence the real function 
R-^ (^^(,2) —^3^(2)) does not have jump discontinuity on [Af^,A^]. Therefore for 2; G M, 
this real function can only have jump discontinuity at the point x*. 

From the coefficient of in fl3.27l) . we see that 

2Re (ef (.)) +if{z) = -^- (a2 + ^) . (3.36) 

(Jj \ iC I 
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Taking the derivative with respect to z, we obtain 



This imphes 

From f l3.27p . it is easy to see that 5^ > as cat < 1. Hence if the derivative of 
{^f (-2) ~ 'C-r(-2)) is non-positive at a point zq G [Af^, A^], then we will have > 0. 

Since Cr{z) is real, this would imply the derivative of the function z defined by fl3.26p 
is positive at the real point m = ^r{zq). By Lemma [II the point zq = z{m) cannot belong 
to Supp (^^N^ = [^i , A^]. This leads to a contradiction and hence we must have 

^Re(ef (^) -e^W) > 0, ze [Af ,A^]. (3.37) 

In particular, if the function Re [^j^{z) — ^^{z)) has more than one zero inside of [A^, A^], 
then at one of the zeros, the derivative in fl3.37p must be smaller than or equal to zero. 
This is a contradiction and hence the function Re (^]^ (z) — ^Ri^)) can vanish at most 
once inside [Xf , A^]. This means that the function Re {^2i^) ~ ^i^i^)) ^ill ^Iso vanish at 
most once inside [A^, A^]. Since we have already shown that Re (^^(-z) — (.^iz)) vanishes 
at least once inside this interval, it must then vanish exactly once inside [A^,A^]. This 
concludes the proof of the lemma. □ 

We can now determine the number of intersection points between and M. 

Lemma 9. The set Sj intersects M at most twice. 

Proof. By Lemma[Hl there exists exactly one point l G [X^ , A^] such that Re (^^(i) — (.si^)^ 
0. Let us assume that Sj and M intersects at a point sq < l. The case when Sq > l can be 
treated similarly. Let us also choose C such that C intersects M at x* < sq. As both Sj and l 
are independent on the choice of C, a different choice will not affect the intersection between 
S) and M. Since x* < Sq, the function Re — is continuous and has different signs on 
(so, l) and (t, 00). However, within these two intervals, the sign of Re — ^|^) remains un- 
changed. Hence there is at most one point z G (sq, 00) such that Re ^ J^^ ^2 ~ ^s^x^ = 0, 

where the integration path is taken along M. Since Re ^ J^m ^2 ~ ^^dx^ = and that S) is 

symmetric with respect to the real axis, we see that there is at most one point on (sq, 00) 
that belongs to Sj. 

Let us now consider the possible intersection points on (— 00, sq). As the choice of 
the branch cut C does not affect the set Sj and its intersection with the real axis, let 
us choose C so that it intersects M at a point Xq > Sq instead. Then Re — ^|^) is 
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continuous on (— cx),So). Moreover, it does not change sign in (—00, sq). Therefore the 

function Re ^ J^n ~ Csdx^ does not vanish in (—00, Sq) ^-nd hence does not intersect 

(—00, So)- This shows that if intersects M at a point Sq < i, then there can at most be 2 
intersection points between S) and M. By using similar argument, one can show the same 
for the case when sq > 

Let us now consider the case when sq = l. If Sq = i, then by choosing C such that 
X* < L {x* > l), we see that Re (^^ — C^^^ does not change sign in (sq, 00) {{—00, Sq)) and 

hence Re (^J^n ^2 " Ca^'^^j does not vanish in either of these intervals. Therefore can 
only intersect M at the point Sq — l. In any case, the set Sj can intersect M at 2 points at 
most. □ 

We can now determine the shape of the set S). 

Proposition 2. The set Sj consists of 4 simple curves, i^^, S^l and S)r. The curve S)'^ 
(^00) open smooth curve that go from (X^ ) to infinity. They approach infinity 

in a direction parallel to the imaginary axis and do not intersect the real axis. The curves 
S)l and Sjji are simple curves joining X^ and X^ . The curve Sjl is in the left hand side 
of S^R in the complex plane and each of these curves intersects the real axis once. These 
two curves are smooth except at their intersections with real axis. Let xl and xr he the 
intersection points of Sjl and 9)r with then {xl, xr) fl [A^, A^] ^ 0. 

Proof. Let the sets So+ and be the intersections of with the upper and lower half 
planes respectively. Then these 2 sets are reflections of each other with respect to the real 
axis. Within the set there are 3 curves i^J, S)f and coming out of the point A^. 
Let us show that these curves are smooth except at A^. Suppose there is a point zq on i^t 
that is not smooth. This means that the function 62 — is not conformal at Zq. Since Jo 
is independent on the choice of the branch cut C, by changing the branch cut if necessary, 
we can assume that both 62 (z) and 9^ (z) are analytic at Zq and therefore the derivative 
of 02 — 0^ must be zero at Zq as the function is not conformal at Zq. This would imply 
^2 {^0) = 'C3^(-2o), which is impossible as the only points where this happens are the points 
A^. Therefore the curves i^^, j = 0, 1, 2 are smooth except at the point A^. 

We will now show that the curves S)^ cannot be connected with one another except at 
the point A^. Suppose the curves Sjj' is connected to Sj'^ at a point Zq ^ X^ . Since both 
curves S)j and S)l are smooth, the curve S)j U ij^ forms a close loop in the upper half 
plane. Let V be the region bounded by this close loop. Then by changing the choice of C 
if necessary, we can assume that the functions 62 {z) and 0^{z) are analytic in the interior 
of V. Then the function Re (^62 — 0^^ is a harmonic function in the interior of V and has 
constant value at the boundary of V. Therefore, by the maximum modulus principle this 
function must be a constant in V. This is not possible and hence the curves Sj'j cannot be 
connected to each other. 

By inspecting the behavior of — at 2; = 00, we see that one of these curves must 
be an open curve that approaches infinity at a direction parallel to the imaginary axis. We 
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Figure 2: The set for a = 0.9, Pn = 0.7 and cat = 0.4. The branch points are given by 
^ 0.12518, ^ 2.48841, A^ ^ 2.40520 + 3.2516i and Af ^ 2.40520 - 3.2516i. The 
point L in Lemma [H] is given by 6 ~ 0.602. For a < 1, the function Re (02 ~ ^si^)) is 
negative in the open region on the left of Sjl and positive in the region Qji on the right 
hand side of S)r. 



will call this curve and its reflection with respect to the real axis Sj^. Since the other 2 
curves cannot intersect each other, and they cannot go to infinity either, they must end at 
the real axis and be connected to the curves in i^". We will call the curve on the left hand 
side and the one on the right hand side S^~^. These two curves must end at different 
points on the real axis as they cannot intersect. Let us denote the curves 9)l and S^r by 

where S^~[ and 9)]^ are the reflections of i^^ and with respect to the real axis and xl, 
xr are their accumulation points on the real axis. 
Let us now show that 

S)+ =Sot,US)tuS)+. (3.39) 

Suppose there is a point zi G Sj^ that does not belong to any of the curves in the right 
hand side of fl3.39p . Then zi must belong to a curve S)^ G S)^. By changing the definition 
of C again if necessary, we see that must be smooth. This curve cannot end on the real 
axis because by Lemma ^ the set Sj can at most intersect the real axis at 2 points and 
has already intersected the real axis at the 2 points xl and xr in f l3.38p . The curve 9)^ 
cannot approach infinity or intersect any other curves in S}~^ either and therefore it must 
be a close loop in the upper half plane. As S)^ cannot intersect the curves in the right 
hand side of f l3.39p . the point \^ must lie outside of the region V bounded by Sj^. This 
would then imply that the harmonic function Re ~ (^3) is constant inside the region 
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V, which is not possible and hence we have 



^ = i5+ U U fiL U ^^j. 



Finally, if {xL,Xf{) fl [\^ , A^] = 0, then the function Re (6*^ — 9^) is harmonic inside the 
region bounded by and S)r, which is not possible as it would imply that it is a constant 
function in this region. This concludes the proof of the proposition. □ 

The shape of the set 9) is indicated in Figure [21 The Octave generated figure shows the 
set for a = 0.9, j3]\f = 0.7 and cat = 0.4. 



3.2.2 Jump discontinuities of the functions 

From now on, we will choose the branch cut C to be a simple curve joining and A^ that 
is symmetric with respect to the real axis. We also require C to lie between the curves l 
and Sjf( in Proposition [2] and that it intersects M at a point A^ < x* < A^. The integration 
contours for the functions 9^{z) in fl3.29p are chosen such that they do not intersect the 
set {—CO, A^) U C and the point Xf in f l3.29p is chosen to be A^. 



Proposition 3. If a > 1, then A^ is a branch point of ^^{z) and is a branch point 
0/^^(2;). On the other hand, if a < 1, then A{^ is a branch point of C,^{z) while X2 
branch point 0/^^(2;). 



^ IS a 



Proof. Let /i = 2, Z2 = 3 for a > 1 and /i = 3, /2 = 2 for a < 1. Then for large enough 
Zq > X2 , all the functions ^j^ are real and we have ^^^(2:0) > ^j^ (zq) > ^j^ (zq) by (11.171) . This 
ordering must preserve at A^ as the roots cannot coincide between zq and A^. At A^, one of 
the roots must coincide with and from the ordering ^i{zq) > (zq) > (zq), we must 
have ^i^(A^) = (A^) = 7^. On the other hand, for small enough < e < X^ , all three 
^j^ will be real. From the asymptotic behavior of C,i{z) at (13.251) . we see that for small 
enough e, ^i{e) will be smaller than both ^|^(e) and ^^{e). From the asymptotic behavior 
of C,2 cind at z = —00 (I1.17P and the fact that these 2 functions has no singularity and 
cannot coincide in (—00, X^), we see that at 2; = e, we must have (e) > ^/^(e). Therefore 
we have (e) > ^/^(e) > ^f'(e). This ordering must again be preserved at X^ . Therefore 
we must have C,i{Xi) = ^/^(Af^) = 7^^. This shows that A^ is a branch point of C,i^ while 
Xi is a branch point of □ 

We can now determine the jump discontinuities of on the branch cuts 



(3.40) 



where A;2 = 1, ^3 = 2 for a > 1 and A;2 = 2, /^a = 1 for a < 1 and Q3j are defined by 

®i = [Af,x*), *B2 = (x*,A^]. (3.41) 
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Figure 3: The branch cut structure of the Riemann surface Cn when a > 1. 



The branch cut structure of the Riemann surface /^at is indicated in Figure [3l 



Let us define 0^ [z) to be constant shifts of the 9^ [ 



z . 



'''^'>-\e^{z)-e^^4X^), a<l. (3.42) 

^^3"^W-^^3^+(Af)+<4Af), a<l. 

Note that the difference between 62 {z) and 0^{z) are the same as the difference between 
e^iz) and e^{z). 

Lemma 10. Let Of {z) he defined as m then we have 9^ (z)- 6^ (z) = e^{z)-9^{z). 

If z E (— cxD, X2) U C, then we also have 

. i (3.43) 

oU^) - eU^) - (eU^) - el^{z)) = o 

where the '+ ' and '— ' subscripts indicates the boundary values on the left and right hand 
sides of (—00, A^) U C. 

Proof. Let li = 2, I2 = 3 when a > 1 and li = 3, I2 = 2 when a < 1. Then from fl3.42p . we 
have 

9^{z) - ~9Uz) - [eUz) - e^z)) = (-1)'^+^ (C+(Af ) - C+(Ar) - C(Af )) • (3-44) 

From (13.421) . it is clear that if 2; G (— 00, A^) U C, then the differences in (13.431) are also 
given by the constant on the right hand side of (I3.44p . 

Let us consider the difference ^i^,+ (Af )-^,^_+(A^). We will compute ^/^',+ (Af ) (^/^,+ (A^)) 
with an integration path that consists of 2 parts. The first part goes along the left (right) 
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hand side of C from to x*. The second part goes along the positive side of the real axis 
from X* to A^ (-^2')- We then have 



C+(Af) 



{x)dx + 
{x)dx + 



Af 



[x)dx, 
{x)dx 



(3.45) 



From (Km . we see that = Cz^,- along C and 



C,i^{x) along (x*,Af), while 



^i^+{x) = ^f_{x) along (x*, A^). From this and fl3.45p . we obtain 



Af 



(3.46) 



where the last inequality follows from the fact that ^i{z) is analytic across C and hence 
we can choose the integration path for Oi{X^) to be along the real axis. 



From dSaSD, we see that ^/;,+(Af ) -^/^,+ (A^) -^?;_(Af ) = and hence 9. 



e. 



N t 



6. 



□ 



From the behavior of ^j{z) on the cuts, we have the following analyticity properties of 
the d^{z). 

Lemma 11. The function 9^ {z) is analytic on <C\{— 00^X2]. The function 62 {z) (9^{z))is 
analytic onC \ ((— c>o, x*] U C) when a > 1 (a < I) and it is analytic on C\ ((—00, A^] U C) 
when a < 1 (a > 1). Let k2 = I, = 2 for a > 1 and k2 = 2, k^ = 1 for a < 1. Let C± 
are the intersections of C with the upper/lower half planes. Then the integrals 0^ {z) have 
the following jump discontinuities. 







= 


+ v±,, ze^k,, j = 2,3, 








= eUz) 


, z e C+, 








= oU') 


+ 2civ/37v7ri, a > 1, z E C^, 






» 


= CW 


— 2cAr(l — /3Ar)7rz, a < 1, z 








= ~el_{z) 


-2civvri, 2G(0,Af], 






[z) 


= oU^) 


— 27ri, z G (—00, 0), 






[z) 


= <-W 


+ 2civ(l — /5Ar)7ri, z G (— 00, 






[z) 


= oU') 


+ 2civ(l - /37v)7rii, z G QSi, 


a < 1 






= CW 


+ 2cNl3NT^i-, z G (-00, Af ], 








= oU') 


+ 2c7v/5Ar7ri, 2; G ^1, a>l 





(3.47) 
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where v± is the constant 

v± = ( K ) - ( K ) + ( K ) - ^1-- ( K ) • (3-48) 



and QSj is defined in ^^3^JJ^. In particular, by the jump discontinuities of the 9^ {z) at the 



points \^ , and the fact that 0^ {z) and Of{z) differs by a constant shift only, we see that 
the constant v± is either —2ciyPiy7ii, —2cn (1 — Pn) m or 0. 

Proof. From the jump discontinuities of ^f{z) in fl3.40p . we see that 

zeQSfe,, J =2, 3 





zeC. 



z e 

The path of integration in the first equation is taken along the real axis and the integration 
paths in the last 2 equations are taken along C. By comparing this with (13.291) and (I3.42p 
and making use of Lemma [TUl we obtain the first two equations in (13.471) . together with 

02,M = OU^) + ^2''(Af ) - 0^{\^.l z e C_. (3.49) 

Let us now compute the constant Q2 (A^) — 6^{Xf). Let li = 2, I2 = 3 when a > 1 and 
h = 3, I2 = 2 when a < 1. Then C,^ will be analytic on (x*, A^] and hence we can compute 
6'/^(A^) using a contour that goes along the right hand side of C. That is, we have 

lN/\N\ I 



erAK)= ei-{x)dx. (3.50) 



where the integration is performed along C. To compute 0^{X^), let us choose an inte- 
gration contour as follows. The integration contour consists of three parts. The first part 
goes from A^ to x* on the left hand side of C. The second part is a closed loop S that goes 
from X* to X* with the branch cut ^2 U C of C,i^ inside it. The last part goes from x* to 
along the left hand side of C. That is, we have 



A4 



C (^f ) = j^^ €4^)dx + iU^)dx + / iZAx)dx. 
Prom the jump discontinuities of if {x) on C, we obtain 
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Figure 4: The contours r+ and r_. 

Since branch cuts of are inside the loop S, this loop can be deformed into a loop around 
z = oo. By computing the integral using residue theorem, we obtain 

02(^4) - ^3^(Af ) = -2mcN{l - (3^), a < 1, 
02(^4) - ^3^(Af ) = 2mcNl3N, a > 1. 

This, together with f l3.49p gives the third and fourth equations in fl3.47p . 
Let us now show that 

e^^^{z) = e^^_{z) + 2cn{i- M^^h ze{-oo,x^]. (3.51) 

The corresponding equation for 9i{z) and 9^ (z) can be proven in a similar way. Let 
z G {—oo, Xi] and let r± be contours from to z indicated as in Figure HI Then we have 

OU^) = [ iUx)dx. 

Let r be the close loop on C such that F = r_|_ — r_, then we have 

= C W + fj^ix)dx, z E (-00, Af]. 

By Cauchy's theorem, we can deform the loop F such that F becomes a loop around ^ = 00. 
By computing the residue, we obtain 

e^^_^{z) = e^^_{z) + 2cjv(l - MTTi, z G (-00, Af ]. (3.52) 

Since 6*^ (z) is a constant shift of 6'^ (z), fl3.52p implies (13.51 1) . 

By using similar argument, we can obtain the rest of the jump discontinuities. □ 

We will conclude this section with the following results on the relative sizes of the 
Re (^9j^(z)^, which are essential in the implementation of the Riemann-Hilbert method. 
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Lemma 12. The real parts of 9^ {z) in jS.^^ are continuous in M+ \ [Ai, A2] and we have 
the followings. 

1. Let xl < xr he the points where intersects the real axis, then we have 

Re - ~e^{z)) > 0, ^ G M+ \ [XL, xrI j = 1, 2, 

Re (e^iz) - 9f^{z)) > 0, ze (0, Af ), (3.53) 
Re{9^{z)-9fl{z))>0, ze{X^,oo). 
where li = 2, I2 = 3 for a > 1 and li = 3, I2 = 2 for a < 1. 

2. Let X* be the intersection point between C and M. Then in a neighborhood D^* of x* , 
we have 

Re {e^{z) - < 0, zeD^*r} C. (3.54) 

where the '+ ' and ' subscripts denote the boundary values at the left and right hand 
sides ofC. 

3. On C, we have 

Re - < 0, zee. (3.55) 

where the '+ ' and '— ' subscripts denote the boundary values at the left and right hand 
sides ofC. 

Proof. From fl3.47p . we see that the real parts of 6f{z) are continuous in R_|_ \ [A^,A^]. 
Now from the proof of Proposition [31 we have, for j = 2, 3, 

Re(ef(^)-^f(^)) >0, ^G(A^,oo), 

Re(ef(^)-ef(^)) <0, ^G(0,Af). 

From the definitions of 9j^{z) (13.291) and 9j^{z) (13.421) . we obtain the second and the third 
equations in (I3.53p . To prove the first equation in (I3.53p . note that if z G {xr, 00), then z 
is in the region Qr in Figure [21 Similarly, if 2; G (O,^^,), then z G Ql- By considering the 
behavior of 62 (z) and 0^{z) at z = 00, we see that 

Re{e!'^iz)-ef^iz))>0, zG(x^,oo), 
Re(^,^(^)-^^,^(^)) >0, zG(0,x^). 

From Proposition [21 we have {xl,xr) fl [Af^, A^] 7^ 0. This implies xl < x* < A^ and 
xr> X* > . Therefore, by the second, third equations in (I3.53p . Lemma [TUl (13.570 and 
the fact that ^i{z) — C,i^{z) is purely imaginary on *Bj, we obtain 

Re (e^iz) - e^^iz)) > Re (9^{z) - 9^^ (z)) > 0, zE (xr, 00), 

Re (e^z) - C (^)) > Re (e^iz) - 9f;{z)) > 0, ze (0, xl). 
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This, together with the second and the third equations in (13.531) . give the first equation in 
(13331). 

We will now prove (13.541) . Let us consider the left hand sides of (I3.54p at the point x* . 
Since {z) — (.[^(z) is purely imaginary on Q3j, we have 

£^ Re (ef (x) - (x)) c/x = 0. (3.58) 

where the integration is performed along the real axis. Therefore we have 

Re (<+(x*)) = Re {el_ix*)) = Re (o^ix*)) (3.59) 

where the '+' and '— ' subscripts indicate the boundary values on the left and right hand 
sides of C. Now note that 6f^j^{x*) is evaluated in the region Vti of Figure [2] and Of^ _{x*) 
is evaluated in the region Again, by considering the behavior of 6*^ (z) and 9^{z) at 
z = oo and using the fact that Re ~ can only change signs across the sets and 
C, we obtain 

±Re (^~^,i(x*) - e^^ix*)) > 0. (3.60) 

Then from (13.591) and (I3.60p . we obtain (13.541) at z = x*. The statement for z G D^.* fl C 
now follows from the continuities of the functions in (13.541) along C. 

Finally, note that inside the set ^2 between C and S^, (See Figure [2]) we have the 
following inequalities. 

Re (ef'^iz) - ^,^(z)) > 0, ze Q2, (3.61) 

by considering z E Q2 on the right hand side of C in (I3.6ip and making use of the jump 
conditions between 62 and 9^ in (13.471) . we obtain (13.551) . □ 

The final result in this section deals with the behavior of these real parts in a neigh- 
borhood of the interval [X^ , A^]. 

Lemma 13. The open intervals {X^,x*) and (x*, A^) each has a neighborhood Ui and U2 
in the complex plane such that 

Re (ef{z) - e^{z)^ > 0, ze Uk^, j = 2, 3. (3.62) 
where k2 = 1, = 2 for a > 1 and k2 = 2, = 1 for a < 1 and 0^ (z) are defined in 

( TO - 
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Proof. Since C,i{z) — ^j^{z) is purely imaginary in ^kj, we have 

/ Re(ef(x)-e,(x))rfx = 0, ze^k^, J = 1,2, (3.53) 

where 03^^ is defined in f l3.4ip . 

Let 05°^, = QSfc^ \ {^kj}- Then on the positive and negative sides of the derivatives 
of the functions 9^.^{z) — Of^±{z) are given by ii^±{z) — ^j^±{z) and are purely imaginary. 
In fact, since ^i,+ {z) = mp^, we see that Cvf (-2^) ~ ~ 2nipN{z) where Pn{z) > 

is the density function of Fn- On the other hand, by the jump discontinuities (13.401) . we 
see that ^i_{z) — C,j^-{z) = —2nip^{z). Hence by the Cauchy Riemann equation, the real 

part of 61 (z) — 0^{z) is decreasing as we move from 03°^. into the upper half plane. From 
(13.631) and (I3.42p . we see that Re (^9^{z) — 9j^(z)^ < for z in the upper half plane near 

'^tj- Similarly, we also have Re (^9^{z) — 9^{z)^ < for z in the lower half plane near 
. This implies (13.621) is true in a neighborhood Uk^ of 5B°_^.. □ 



4 Riemann-Hilbert analysis 

We can now implement the Riemann-Hilbert method to obtain the strong asymptotics for 
the multiple Laguerre polynomials introduced in Section [2] and use it to prove Theorem [3l 
The analysis is very similar to those in [T2] (See also [26]). 

Let C{f) be the Cauchy transform of the function f{z) G L^(R4_) in M_|_ 

CifKz) = ^ [ ^ds, (4.1) 



27ri 

and let Wi{z) and W2{z) be the weights of the multiple Laguerre polynomials. 

Denote by ki and k,2 the constants 

Ki = -27Ti f /i^-i /V, ) , 1^2 = -2m (hf^ 



''Nq-1,NiJ ' "-2— y''NQ,Ni-l 

Then due to the orthogonality condition (12. ip . the following matrix 

/ Pno,Ni{z) C{Pno,N^Wi){z) C{Pno,NiW2){z) 

Y{z) = KiPno-1,Ni{z) KiC{Pno-1,NiWi){z) KiC{PNa-l,NiW2){z) | (4.3) 
\I^2Pno,Ni-i{z) K,2C{Pno,Ni-1Wi){z) K2C{Pno,Ni-1W2){z) ^ 

29 



is the unique solution of the following Riemann-Hilbert problem. 

1. Y{z) is analytic in C \ IR+, 

/I wi{z) W2iz)\ 

2. Y+{z) = Y_{z) 1 , z G R 

Vo 1 ; 

;^ 

3. Y{z) = {l + 0{z-^)) I z-^^ 



(4.4) 



z ^ oo. 



z-^' 



4. Y{z)=0{l), z^O. 



By a similar computation as the one in [TU] and [Hj, we can express the kernel (12. 5p in 
terms of the solution of the Riemann-Hilbert problem Y{z). 

{xy)"^ [e-'^y [Y^\y)Y^{x)\^^ + e~''^-'y [Y^\y)Y^{x)\ .\ 

KM,N{x,y) = — — 7 , 

27rt[x — y) 

_ 1\ (4-5) 

(xy) ' 



e-^'y e-^'''''y^ Y-\y)Y4x) lo 



2m{x — y) 

where ^421 and A^i are the 21^^ and 31*'' entries of A. 



4.1 First transformation of the Riemann-Hilbert problem 

We should now use the functions 0^ {z) in (13.421) to deform the Riemann-Hilbert problem 
(14. 4p . Our goal is to deform the Riemann-Hilbert problem so that it can be approximated 
by a Riemann-Hilbert problem that is explicitly solvable. As in [12] and [26], the set S) in 
(13.311) will be important to the construction. 

We will now start deforming the Riemann-Hilbert problem (14. 4p . First let us define the 
functions gf{z) to be 



g^iz) = e^iz) + (1 - cj,) \ogz, g^{z) = 6^ (z) + z, 



z 



(Jj 



(4.6) 



where 0^ (z) is defined in (I3.42p and the branch cut of log 2; in gi{z) is chosen to be the 
negative real axis. 

We then define T{z) to be 

T{z) = dtag (e-*^'", e"^^'", e"*^'") Yiz)dtag (e*^<(^), e^'^^^^\ e^'^^^-'^) , (4.7) 
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where I2 and are given by 



JN 

to 



^2^-^2'f+(Af)+^^(Af), a>l; 

rN_S a>l- 

il!{z)-ei'^4x^) + e^^_{x^), a<i. 



The matrix T{z) will satisfy the following Riemann-Hilbert problem. 

1. T{z) is analytic in C \ (M U C), 

2. T+(2) = T_(2)Jt(z), ^gMUC, 

3. T{z) = I + 0{z-^), z^oo, 

4. T(z) = 0(l), z^O. 

where Jt{z) is the following matrix 

Mz) 



(4i 



^gM(ef^^{^)-ef^_(^)) ^M{e^^^(z)-e^_(z)) ^M{ei'_^iz)-e^_iz))\ 
\ e*^(^'3:+W-e'3:-{-))y 



(4.9) 



By applying Lemma fTTl to the 9^{z), we can simplify the jump matrix Jt{z). In particular, 
on ®fe2 ill f!3.4ip . we have 



f^M{e^(z)-~9^(z))_ 






1 e^^(^'3:+{^)-eT-{^))\ 
M{e«{z)-e«{z))_ Q 

1 / 



(4.10) 



while on 23 ^3, we have 



Mz) 



/ M{e^(z)^~e^(z))^ M{e^^4z)-el_{z)) 



\ 








1 





1 



M{e^{z)-e^{z))_ 



(4.11) 



On the rest of the positive real axis, the jump matrix becomes 
Jt{z) = 1 

Vo 1 



(4.12) 



This is because McnPn = Ni, Mcn{1 — Pn) = Nq and Mcn = N are all integers. And 
on the negative real axis, the matrix T{z) has no jump for the same reason. Note that the 
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Figure 5: The contours El and Er. 



jump matrix Jt{z) is continuous at 2; = as the off-diagonal entries of fl4.12p contain the 
factor which vanishes at the origin. 

The jump matrix on C is given by 



/I Q \ 

^0 ^ 



(4.13) 



The Riemann-Hilbert problem for T{z) now takes same form as the one in pjj (See also 
[26] ) and the techniques developed there can now be applied to our problem. 

4.2 Lens opening and approximation of the Riemann-Hilbert 
problem 

We will now apply the global lens opening technique developed in We will have to 
use the properties of the set ^ in (13.311) to define the global lens contour. 

By Proposition [21 the set union of the set and C divides the complex plane into 4 
regions VLr-, ^1 and ^2- The region VLl {^r) is an open region that lies on the left 
(right) hand side of the contours S)^ and S)l. {S^b) The region f^i ((^2) is the region 
bounded by S)l {^r) and C. (See Figure [2]) 

By considering the behavior of 62 (z) — 0^ {z) near 2; = 00 in (13.301) . we see that for 
a < 1 (a > 1), Re {O2 {z) — 0^{z)) is negative (positive) in VIl and positive (negative) in 
Qji- Let us define the contour (S^) to be a contour from to in (^i?) such 
that El intersects M at a point xi < and E^ intersects M at a point X2 > A^. (See Figure 

E> 

Let us define the lens contours H^, j = 1, 2 around the branch cut [X^ , A^] as follows. 
The contours S^. are contours in the neighborhoods Uj in Lemma [13] joining Aj^ and C in 
the upper/lower half plane. Together with the contours El and Er, the lens contours are 
depicted in Figure [H] Let the matrices Qi, Q2, Q3, Qi be the followings 
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Figure 6: The lens contours. 

'10 \ /I 0^ 

Gi={0 1 _eA<^l(-)-^T(-)) ], g,= \0 1 

.0 1 VO _e^<^TW-^?W) 1^ 

(4.14) 

1 0\ / I 0^ ^ ^ 

^3 = ( _^M{e-i.)-e-i.)) 1 , ^4 = 10 

1/ \_^M{e-i.)-e-i.)) Q 

We will now define the matrix S{z) to be S{z) = T{z) outside of the lens region. Inside 
the lens region, we define S{z) to be 



(4.15) 



For a < 1, S{z) = T{z)g2, z G Di, S{z) = T{z)gi, z G D2, 

S{z)=T{z)G2Qt\ zeLl, S{z)=T{z)G^Gf\ zeLl. 
Fora>l, S{z)=T{z)gi, zeBi, S{z)=T{z)g2, 2 G D2, 

siz) = T{z)g,gt\ zeLl, siz) = T{z)g2gt\ zgl|. 

Let S be the union of the lens contours. 

E = ELUEnU {UU^I U Ej) . 

Then by using Lemma [TTl it is easy to check that S{z) satisfies the following Riemann- 
Hilbert problem. 

1. 5(;z) is analytic in C\(M+U sue), 

2. S^{z) = S4z)Js{z), ^G(M+UHUC), 

3. S{z) = I + 0{z'^), z-^00, ^ ' ' 

4. S{z) = 0{l), z^O. 

where the matrix Js{z) is given by the following. Let C = U'j^qCj where Co is the boundary 
between Di and ©2) Ci is the boundary between and and C2 is the boundary between 
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L\_ and L^. Then on C, the jump matrix Js{z) is given by 

AO \ 

For a < 1, Jsiz) =00 1 , z G Co, 

\o -1 e^^(^'3':+w-«i-w)y 

/I 

= 1 

zeCi, 1 = 1,2. 

n o\ 

For a > 1, Jsiz) = Q e*«+(^)-'^(^)) -1 , z e Cq, 

\0 1 0/ 

/I 

\ 10 

zeCi, 1 = 1,2. 

On the lens contour S, the jump 75(2;) is given by 



For 


a < 1, 


Jsi 




= 02 


1 


z G 




J5( 




= 6^1, 


^ G Er, 






Jsi 




= ^4- 


1 


z G 


"1 ' 




[^) 


= ^3-\ 


zeEf 


For 


a > 1, 


Jsi 




= ^r 


1 


z G 


Si, 


Jsi 




= G2, 


z G Sij, 






Jsi 




= ^3 


1 

5 


z G 


"1 5 


J5I 


[^) 


= 


zeEf 



(4.17) 



(4.18) 



Finally, let Mj be the following. 

Mi = (0,Af), R2 = iX^,X2). (4.19) 
Then on the real axis, the jumps of S{z) are given by 

A e*<^'2':+(-)-ef:-(-)) o\ 

^5(2) = 1 , 2 G Rk2, 

\0 1/ 

/l e*^(^'3'^+(")-^'^-("))\ 
^5(^2) =01 , ^ G Mfc3, (4.20) 

\0 1 / 

1 0\ / 1\ 

Jsiz) =1-10 0, ze Jsiz) = 10, zE Q3fc3. 

1/ \-l 0/ 

where k2 = 1, = 2 for a > 1 and /c2 = 2, fcs = 1 for a < 1. On the rest of the positive 
real axis, the jump is given by fl4.12p . 
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By considering the behavior of O2 {z) — 9^ {z) near z = 00 in (13.301) and by making use 
of Lemma [TUl we see that for a < 1 (a > 1), Re {o^ {^) ~ ^si^)^ is negative (positive) 
in Si and positive (negative) in S/j. Hence, from (14.181) . we see that the off-diagonal 
entries in Js{z) become exponentially small on and as M — > 00 outside of some 
small neighborhoods Dk around the points A^. By Lemma fT3l we see that, away from the 
points A^, the jump matrices on the lens contours Ef and Ef are all exponentially close 
to the identity matrix as M ^ 00. By (I3.53P and the fact that X2 > xr, it follows that 
Jsiz) are also exponentially close to the identity matrix on M+ \ ([A^, A^] U Di U D2) as 
M ^ 00. Now by (ESD, we see that the for a < 1, the 31*'' entry of Jsiz) in C1UC2 ( 14T71) is 
exponentially small in M, while for a > 1, the 21*^ entry of Js{z) in C1UC2 is exponentially 
small in M. Finally, by (I3.55p . we see that for a < 1, the 33*'^ entry of Js{z) in C is also 
exponentially small in M, while for a > 1, the 22*^* entry of Js{z) in C is exponentially 
small in M. This suggests the following approximation to the Riemann-Hilbert problem 

dHSD. 

1. S^{z) is analytic in C \ ([Af , A^] U C) , 

2. S^iz)=S^iz)J^iz), zG[Af,A^]UC, 

3. S°°iz) = I + Oiz-^), z^oo, (^-21) 

4. 5°°(^) = O ((^ - Af )-i) , z^Xf, j = l,...,4. 

where the matrix Joo{z) is the same as Js{z) on [Af , Af ] and on C, it is given by 

'1 0^ 



For a < 1, Jooiz) 



For a > 1, Jao(z) 




(4.22) 



In the next section we will give an explicit solution to this Riemann-Hilbert problem and 
we will eventually show that S^{z) is a good approximation of S{z) when z is outside of 
the small neighborhood Dk of the branch point Af . 



4.3 Outer parametrix 

Let be the Riemann surface defined by (13.271) and let Tj be the images of QSj on 
under the map ^f^(z). That is 

r, = {(z,Oe£^| e = ef+W, ze^,.}, j = i,2. (4.23) 

Similarly, we define Fc to be the image of C. 

r,= {(^,Oe£^| e = C+W, ZEC.}, (4.24) 
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where h = 2, I2 = 3 when a > 1 and h = 3, I2 = 2 when a < 1. 

Let us now define the functions S^{^), /c = 1, 2, 3 to be the following functions on ^at. 



sr({) = «, 



\ 1=1 



aVnU(l + 7f) 



(e + i)(£+0 

n,ti(e-7f)' 



a — 1 



n;=i(i+a7i 



(4.25) 



nL(e-7^ 



where 7^ are the roots of polynomial f l3.28p . The branch cuts of the square root in fl4.25p 
are chosen to be the contours Fi, r2 fl4.23p and Tc fl4.24p in 

By using the asymptotic behavior of the functions ^^(-2) (H-l^p and 03.131) . with c and 
/? replaced by and P^, we see that the all the functions 5*^(0 are holomorphic near 
^^(0) for m = 1,2,3. Moreover, at the points ^^(00), these functions satisfy 



SrieZi^)) = Smk, k,m = 1,2,3. 
Let us define S°°{z) to be the following matrix on 2; G C 



(4.26) 



(4.27) 



Then, since the branch cut of the square root in fl4.25p are chosen to be Tj and Tc, we 
have, from the jump discontinuities of the f l3.40p . the following 



= 1,2,3, / 
1,2,3. 



2,3, 



(4.28) 



From this and the asymptotic behavior f l4.26p of the S'?°(^) at z = 00 and its behavior at 
the branch points Aj^, we see that the matrix S°^{z) satisfies the Riemann-Hilbert problem 

dEni). 

Proposition 4. The matrix S°°{z) defined by ^4-21 ) satisfies the Riemann-Hilbert problem 



4.4 Local parametrices near the edge points 

Near the edge points A^, the approximation of S{z) by S°°{z) failed and we must solve the 
Riemann-Hilbert problem exactly near these points and match the solutions to the outer 
parametrix (14.271) up to an error term of order 0{M^^). To be precise, let 5 > and let 
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Figure 7: The contour T^. 



Dk be a disc of radius 6 centered at the point A^, k = 1, . . . ,A. We would hke to construct 
local parametrices S'^(z) in such that 

1. S''{z) is analytic in \ (D^ n (M U H U C)), 

2. S^{z) = St{z)Js{z), z e Dkn{RUEUC) , (4.29) 

3. S^z) = (/ + 0{M-^)) S°^{z), z G dDk. 

The local parametrices S^{z) can be constructed by using the Airy function as in [12] (See 
also [26]). Since the construction is identical to that in [12] and [26], we shall not repeat 
the details here and refer the readers to these 2 papers. 



4.5 Last transformation of the Riemann-Hilbert problem 

Let us now show that the parametrices constructed in Section 14.31 and Section 14.41 are 
indeed good approximation to the solution S{z) of the Riemann-Hilbert problem f l4.16p . 
Let us define R{z) to be the following matrix 

^ I S{z) {S''{z))'\ 2 inside Dfc, A; = 1,...,4; . 
^' \ Siz){S^iz)y\ 2 outside of Dfc, A; = 1,...,4. ^' ' 

Then the function R{z) has jump discontinuities on the contour Tji shown in Figure [71 In 
particular, R{z) satisfies the Riemann-Hilbert problem 

1. R{z) is analytic in C \ Fij 

2. R+{z) = R4z)Jr{z) 

4 31) 

3. R{z) = I + 0{z-^), z^oo, ^ ' ' 

4. R{z) = 0{l), z^O. 

From the definition of R{z) (14.301) . it is easy to see that the jumps Jr{z) has the following 
order of magnitude. 

J (,) = S ^ + 0(M-i), zedD,,k = l,...A; .4 32^ 

' \ I + O (e"*^'') , for some fixed r] > on the rest of Tr. ^ ' ' 
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Then by the standard theory, [T7], [H], [19], we have 



uniformly in C 

In particular, the solution S{z) of the Riemann-Hilbert problem (I4.16p can be approx- 
imated by S°°{z) and S''{z) as 

^..^(ii + OiM-'))S''iz), zeD,,k = i,...A; 

\ (/ + 0(M-i))S°°(2), z outside of Dfc, A; = 1,...,4. ^ ' 



5 Universality of the correlation kernel 

The the universality result Theorem [3] can now be proven by using the asymptotics of the 
multiple Laguerre polynomials obtained in the last section. The proof is the same as the 
ones in Section 9 of [11], [12] and [28] and we shall leave the readers to verify the details. 



6 Conclusions 

This paper complements the results obtained in our earlier paper ^28j on the universality of 
complex Wishart ensembles 2 with distinct eigenvalues and the number of each eigenvalue 
becomes large as the size of the ensemble goes to infinity. By using results in the Stieltjes 
transform of the limiting eigenvalue distribution, we were able to overcome difficulties in 
determining the sheet structure of the Riemann surface (11.51) . This, together with an 
analysis of the topology of the zero set of a function h{x) fll.l8p . enables us to apply 
the Deift-Zhou steepest descent method to obtain asymptotic formula for the eigenvalue 
correlation function. Together with [28], we showed that unless the discriminant of the 
polynomial (11. 6p is zero, the local eigenvalue statistics are given by the sine-kernel (I1.12p 
in the bulk and the Airy kernel (ll.lSp in the edge of the spectrum. We have also shown that 
the largest eigenvalue is distributed according to the Tracy- Widom distribution (11.160 . 

When the discriminant A of (ll.6p becomes zero, the ensemble goes through a phase 
transition in which the support of the limiting eigenvalue density changes from 1 interval 
to 2 intervals. This phase transition can be studied by the method developed in [13] and 
the eigenvalue correlation function will be given by the Pearcey kernel near the critical 
point of the spectrum, the point where the support is splitting. 

The use of Stieltjes transform to provide a Riemann surface needed for the implemen- 
tation of the Riemann-Hilbert analysis can be generalized to cases where the covariance 
matrix has more than 2 distinct eigenvalues. However, when the Riemann surface has 
more than 2 complex branch points, the determination of the zero set in (I3.3ip is very 
complicated and a full generalization to these cases is still a very challenging open problem. 
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